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FOREWORD 

This work represents an extension of the work reported in Naval Surface 
Warfare Center (NSWC) TR 86-253. The latter report presented computational 
methods for predicting aerodynamic loading for supersonic Mach numbers on non- 
axisymmetric flight vehicles at constant or steady incidence. The extension is for the 
prediction of aerodynamic loading or dynamic derivatives associated with constant 
body axis angular rates and/or acceleration. The resulting computer program allows 
one to predict roll and pitch damping, fin Magnus, and other dynamic derivatives for 
the preliminary and intermediate design stage. 

Support for the work was provided by the following sponsors: 

1. Aerodynamics and Structures Block of the Surface-Launched Weaponry Tech¬ 
nology Program Project Numbers RA11G13/RU11G14. 

2. Air-Launched Antisurface Weaponry Technology Project Number RU11811. 

This report was reviewed and approved by Dr. T. J. Rice, Head, Aeromechanics 
Branch and C. A. Cooper, Head, Missile Systems Division. 
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1.0 INTRODUCTION 

Preliminary design requires the estimate of aerodynamics for a large set of 
free-stream and geometric variations. To keep computational costs reasonable, 
rapid (but reasonably accurate) methods are sought and utilized. 

For cruciform finned, axisymmetric bodies at low angle of attack, component 
buildup methods 1 - 2 can be utilized for a restricted range of configurations and free- 
stream conditions that are applicable to current designs. 

Linear surface singularity methods 3 - 4 have been highly developed for complex 
configurations, including high angle-of-attack vortex modelling. Set-up and run 
times are fairly long. 

The second-order Van Dyke model5, which corrects for compressibility effects, 
was extensively modified and adapted for noncircular bodies with planar discon¬ 
tinuities, including inlets, and multi sets of fins. 6 9 The NANC code 1612 includes 
the second-order potential model plus a local solution model based on the methods of 
Reference 13. 

Dynamic derivative estimates are calculated in the NSWC Aeroprediction code 
of Reference 1. Body-alone contributions are given by an empirical data fit. Fin- 
alone contributions are given by thin-wing theory. Interference modelling is incom¬ 
plete when compared with the static case. Thin-wing theory is based on conical 
solutions that are not valid for small aspect ratio fins. 

Dynamic derivative estimation is primarily of interest for unguided appli¬ 
cations. Pitch damping is usually ignored for guided applications since the autopilot 
provides control surface moments, which are proportional to the pitch rate and 
significantly larger. 

Classical unsteady aerodynamics is based upon the unsteady linear potential 
equation. 14 The unsteady potential equation may be applied to harmonic analysis of 
the rigid motion or to aeroelastic applications. The harmonic gradient method of 
Reference 15 is a singularity collocation code utilizing a complex potential. The total 
pitch damping for low frequency is the same as for the constant pitch and accel¬ 
eration rate estimate. 

The NANC code is modified in this report to compute forces and moments 
which are a function of constant axial rotation rates and moment center accel¬ 
erations. The remainder of the report assumes that the reader is familiar with the 
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earlier work summarized in Reference 10 and concentrates on the dynamic 
derivative modifications. 


2.0 GEOMETRY, KREE-STREAM VELOCITY 
AND FORCE CONVENTIONS 


Assumptions concerning body and fin geometry are the same as in Reference 
10. Figures 1 through 3 are taken from Reference 10. However, the body is assumed 
to have a pointed nose for the low Mach number range. 

A static blunt body model establishes a matching plane of velocity component 
data by matching a modified Newtonian pressure distribution and utilizing conical 
potential functions and other assumptions. A blunt body model for the plunging 
acceleration problem was not deemed feasible. A model for pure axial rotation is 
only somewhat more feasible. 

Computations utilize a cylindrical coordinate system and the "thin-fin” approx¬ 
imation with the fin midplane lying on a cylindrical coordinate ray. The body is 
divided into sections by planar discontinuities (including planar inlets). Any section 
except the first may have a set of fins with up to six fins. See Reference 10 for a more 
detailed discussion. 

Force convention is the same as used in earlier work. The axial force, Fa, acts 
in the x direction; the normal force, F\, acts in the z direction; and the side force, Fy, 
acts in the negative y direction. The roll moment, Mf, acts in the negative x 
direction; the yawing moment, M n , acts in the negative z direction; and the pitching 
moment, M m , acts in the negative y direction. 

The main difference between the current work and earlier work is the body- 
axis oriented equivalent dimensionless free-stream velocity vector. In Cartesian 
coordinates, the velocity vector relative to body axes, in dimensionless form, is given 
as 




V -q'z+r'y 

— + - + cos« cosp 

V V 


+ 


V — r'(x — x') + p'z 

— + - + sinG' j' 

V V 


( 2 - 1 ) 


+ 


w o‘ p v + q (x - x ) 

- + -+ Sina cosp 

v V, 


k' 


‘I, > q,« 


lim t —* 0 
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Note that the time derivatives of the axial rotation rates could be included as well, 
but are neglected, p', q\ r' are the axial rotation rates in the negative x, y, and z 
directions, respectively, and uo, vo, and wo are the corresponding acceleration rates. 
For the total pitch damping problem, q' and wo are the pertinent input terms. « and 
P' are the angle of attack and side slip, respectively. 

An additional symmetry operation variation for LM (defined in Reference 10) is 
introduced here. LM = 1 is for axial force as the only nonzero force term. LM = 2 is 
for axial force and rolling moment as the only nonzero forces. The new symmetry 
value for LM = 2 is for a circular body with quarter-plane symmetric fin distribution 
and a roll control deflection or constant roll rate. For the LM = 2 case, the 
computation is LM = 1 before the first fin set and LM = 2 subsequently. LM = 3 is 
for the axial force, normal force, and pitching moment nonzero. LM = 5 is for axial 
force, side force, and side moment nonzero. LM = 6 is for a full plane computation. 
LM is defined for each section of the body and must increase or remain the same as 
the computation is marched down the body. 


3.0 THEORETICAL DEVELOPMENT 


3.1 FIRST- AND SECOND-ORDER POTENTIAL EQUATIONS 
The full nonlinear, dimensionless, potential equation is 


a~v • q -- \i; 


h J_ Q VQ 2 
V“ rft 2 V , dt 2 


(3-1) 


Here, a is the speed of sound divided by a, and Q = is the dimensionless velocity 
vector. Equation (3-1) is appropriate to a body moving through a fluid at rest. It is 
derived from the continuity equation, momentum equations, and the Bernoulli 
relationship. For body-fixed coordinates, 


<) 

- - V-V 

(3-2) 

d t 

(3-3) 

"V, 


In addition, the relative velocity, q r = Q - V, is introduced. The continuity 
equation then becomes 


— - — + V • q + - (q • V)p = 0 (3-4) 

V r p (It ’ p r 
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The Bernoulli relationship for body coordinates becomes 


f 1 cXj> 

\I“ — — + 
r v at 


q; - v- 


(3-5) 


V - 1 V - 1 

The density may be eliminated from Equation (3-4) by using Equation (3-5) to yield 


a'-V -q = \r 
1 r * 


• V 


\ ~ V " » 1 o 2 4> a / q^ “ v " 


1 a 4> a / M , \ aQ 
7 77 + 7 ( 9 ) + V 77 


(3-6) 


V" at 


v - 1 ■> 

a" = 1 + --- \1- 


cxb 

2 ( —) 


V“ - q‘ 


(3-7) 


A small disturbance from the static free-stream dimensionless velocity vector is 
given by 


t = »’ + (c w - r >+ q = >' + q'*+ q 


(3-8) 


Substitution of Equation (3-8) into Equation (3-6) and neglecting higher order 
terms yields the first-order wave equation 


V-q\=3 


M 


v 2 at 2 


+ 2 


8u i \ + d 2h 

' ax 2 


at 


(3-9) 


q', is the first-order velocity vector, u' the x component of qj, and 4>j the corresponding 
potential function. 


The first-order problem may be further broken down into three steady prob¬ 


lems. 


M; 

♦i = + ~i + v . 


X M; 

p 2 v 


j*b 


(3-10) 


Here 4>i is the potential function associated with the angle-of-attack problem and 
axial rotation rates. <J> a and $b are equivalent steady potential functions associated 
with the plunging rate problem. Equation (3-10) satisfies Equation (3-9) when <J>i, 
<J) a , and 4>b satisfy the steady first-order wave equations 


2 

<V. q ) [T— =0 

f c cTx 


(3-11) 
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(V x q.) = 0 (3-12) 

f c 

( = 1, a, b; c stands for the crossflow plane and components. 

As in earlier work, an improved solution to the first-order problem is obtained 
by evaluating the neglected nonlinear terms using the first-order solution and 
solving a nonhomogeneous potential equation. Only time-independent terms are 
considered. 


( V-q,) c 

- p 2 — = \i 2 v • a 

dx 

(3-13) 

(V x q.,)^ 

= 0 

(3-14) 

il = U 

i’ + il 

(3-15) 


Q — u. 
x i 


1 2 ( 2 - \ ) ,4 2 ~ \ ) <2 2 , 

-q - - M u, — - M q. - 1 

o *1 a l 9 x 


Y — 1 •> 2 

- M' u; + q. +q • q, 

9 » l M l 


(3-16) 


Q c = <»,, 


<*1 2 - 


1 \t'2 ..2 


+ -- M- u; 

2 2 ” 1 


+ q <q7 + q -q,> 


M 


+ — (y — 1) 

r 


V + ( 'k X(q iC Xq j 


(3-17) 


Boundary conditions for first- and second-order problems are 

(q + q )• n = 0 (3-18) 

e r = 1,2 

The equivalent steady problem boundary conditions for tiie q^ problems are 


(vr)- = 


* 


q + x t - 

• n = - — <t> 

a V / 

flX 

' q. 

q ;1 f x ~ ) 

l-n -- -W 


body 


fin 


(3-19) 


(3-20) 
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Here, dr ,/dx and t x are body and fin slopes, respectively. 


3.2 COMPUTATIONAL COORDINATES AND GRIDS 

The region between the body and the Mach cone is mapped to the rectangular 
region shown in Figure 2 by the shearing transformation 



Further clustering transformations i, = HO and <t> = <t>(B) are used for nonuniform 
gridding of the B and £, variables. The functional dependencies are not given 
explicitly. See Section 3.2 in Reference 10 for more details. 


3.3 NUMERICAL METHODS 

Most of the numerical methods are as reported in Reference 10. However, some 
of them have been changed. 

The implicit formulation for the first body section is as in earlier work. The 
velocity vector advancement Equations (3-21) through (3-14) are single second-order 
equations in a potential, <\> = xF.f." F is known as a conical potential function. 
However, actual conical similarity for the total pitch damping problem requires the 
functional form of <£ = x (G + xH). The pressure distribution on a cone is now linear 
with x instead of constant. Therefore, the computation uses more than two steps for a 
solution for a cone. The solution at x = 0 is assumed to be conical. However, the 
resultant computation for a numerical solution using a few marching steps does not 
vary significantly from a solution developed based on the true conical similarity of <J) 
= x (G + xH) and two marching steps. 

At body planar discontinuities and supersonic leading and trailing edges, the 
jump in various velocity components is obtained by application of the method of 
"weak solutions” 1 r > combined with a downstream boundary condition and conser¬ 
vation of the velocity component tangent to the discontinuity edge condition. At a 
subsonic leading edge, the solution for the velocity jumps does not exist and the 
conservation relationship provided by the method of "weak solutions” must be 
replaced by a heuristic one. The "weak solution” conservation relation is modified 
due to the modified Equation (3-13). 

For all sections, except the first, the conservation velocity vector advancement 
equations are solved using a MacCormack 1 ? predictor-corrector scheme for points not 
on a solid surface. Body and fin surface velocities are advanced by characteristic 
compatibility relations (two in number) combined with a solid surface boundary 
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condition.18 At a corner where the fin and body meet, no unique solution is possible. 
The advancement equations for the fin-body junction have been modified. The fin- 
body junction line is assumed to be a streamline. Therefore, fin and body boundary 
conditions and one other condition are needed to solve for the three velocity 
components. 

One estimate of the axial velocity component, u, can be obtained from one body 
advancement quantity and the body boundary condition. 


e b = u bP +t 0 + V B- e 0 W B 


V B ~ C 0 W B = (U B + u ' j - v » + £ 0 W - 

Here, u, v, and w are the axial, radial, and 0 components of velocity, ey = l/rb(drb/<30). 
The B subscript stands for the body. For the first- and second-order problem, ui, v;, 
and w' r are free-stream velocity components. For the <J>b problem, q’ x = q*. For the <J> a 
problem, u' r = u r x/V„ - $b, v; = v*x/V r , and wj = w*xAAo. u B is obtained from Equa¬ 
tions (3-22) and (3-23). E b is known. 

A second estimate of the axial velocity component is obtained from the fin 
advancement quantity and the fin boundary condition 


E f = ± Pu K + w f . 


w t , = - w' + t (u_ + u'J 

p / x F r 


Here, Ef’ is known, up may be obtained from Equations (3-24) and (3-25). The upper 
sign is for the 9 > 0f, n side of the fin. 

The final values of u, v, w for the fin-body junction are obtained from 


u = - ( U B + Up) 


v - cw = — (u + u’) — v'+ e w' y 

dx 0 

w = — w' + t (u + u' r ) (3-28) 

r X 

Note that $b need only be determined on the solid surfaces and is needed for the <J> a 
problem boundary condition and evaluation of the Bernoulli pressure coefficient 
relationship. An advancement equation for <J>b is given by 


— < 1 - £,) + - v + t w 
cTx (1 b x 1 
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On a body surface or interior point, t x is set to zero. 

A subsonic leading edge requires a modification of the jump relations. For all 
cases, 


k.j 




k-I.j 


2 2 
r 0 “ r k.j 


(3-30) 


k is the grid index in the r direction;] is the index in the 0 direction. Equation (3-30) 
combined with downstream boundary condition and conservation of the velocity 
component tangent to the edge provides a solution for the downstream values of 
velocity. The square root ratio is limited to a value of 2. ro is the radius of the 
leading edge. Equation (3-30) has the well known square root singularity for a 
subsonic leading edge. 

3.4 LOADING COEFFICIENTS AND OTHER NUMERICAL CONSIDERATIONS 


The axial acceleration rates contribution to velocities in the limit as t -* 0 is 
given by differentiation of the last two terms of Equation (.3-10). 


M; 

q.„. = —7 


q„ - XC L 


- ■ 7 


(3-31) 


The first- or second-order pressure coefficient is then given by 


C p =(j 3 - 5 - 1|/(.7M;) 

(3-32) 

x= 1 + .2M;(q';-Q^- 2 4> b ) 

(3-33) 

% = T,* + ( U + q <’ 

(3-34) 

€ = 1,2 



This is a "hybrid” model for the plunging rate case since q ax is a first-order potential 
quantity. 

The inviscid loading coefficients are as given in Section 5.0 of Reference 10. 

Smoothing of the MacCormack vector quantities is an input option for sections 
with fins. It is particularly needed for subsonic leading edges. Q|, with q ax = 0 is 
used as a weighting function for first- and second-order vector terms as in Section 
3.11 of Reference 10. Q|, with q ax * 0 is used as a weighting function for the <J) a and 
4>b problem vector terms. 
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The explicit marching solution will fail when X becomes negative in Equation 
(3-32). For q, B = 0,X is set to zero or the vacuum value and the velocities are adjusted 
as discussed in Section 3.11 of Reference 10. For q* * 0 and X negative, the solution 
is halted and an error message is written to the output File. 


3.5 HIGH MACH NUMBER SOLUTION 

The potential model breaks down at Mach numbers where the origin Mach cone 
crosses the body surface, q,. is set to zero for this local solution model. 

The velocity vector, q*, is normalized to 1 as 

q»=q^|q«l (3-35) 

q^ is then used to find the turning angle between an effective free-stream velocity 
vector and the local surface normal. These are combined with the local solution 
methods of Reference 13. It is expected that these local solution methods will be 
inaccurate below M*. = 4. Note that blunt bodies may be considered. 


4.0 EVALUATION OF THE NUMERICAL METHODS 


Evaluation is almost exclusively made by comparison with experimental data. 
The data sources are ballistic range data or wind tunnel tests. Dynamic derivatives 
must be based on parameter estimation techniques associated with kinematic data 
and an assumed aerodynamic model. In general, the accuracy and repeatibility of 
the estimated coefficients is much worse than for static coefficients. 

The code development is for nonaxisymmetric bodies with combined angle of 
attack, control deflection, rotation rates, and axial accelerations. Data, however, is 
mostly available for C| p ' and C m q' + C m i for circular bodies at zero incidence. 
Computations are separate for C[ p ', C mq ', and C m d. 

The dynamic derivatives are given in general form as 


' ^H X R < ^D (a> i X D /V f ) 


(4-1) 


Sr is a reference area; xr is a reference length; M; is a force or moment; Qd is the 
dynamic pressure; xo is a length associated with a reduced frequency; and co; is an 
axial roll rate or dimensionless acceleration. 
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4.1 BODY-ALONE COMPARISONS 

The reference area for bodies will be the maximum cross-sectional area. 

The first computation is for a 2.98 caliber length cone. The moment center is 
located 2.18 calibers from the nose, is the maximum body diameter and x[) is the 
maximum body radius. Unless otherwise indicated, references will not vary in the 
rest of the report. Figure 4 compares C m q' + Cm* computational values with values 
extracted from the ballistic range data of Reference 19. The usual lack of repeat- 
ibility of the data is shown due to different initial yaws and, hence, epicyclic history. 
Figure 5 compares C m q' + C m i computational values for a cone-cylinder and the 
ballistic range data of Reference 19. 

The next comparison computations are for the Army-Navy Spinner config¬ 
urations. These configurations have 2-caliber secant ogive noses with arc radii twice 
that of a tangent ogive and various body lengths and moment center locations. 
Figure 6 compares C m q' + C m i computations with data for a 5-caliber length body. 
Note the great differences between different range tests and wind tunnel tests. Also 
plotted are C m a and a slender body value. Figure 7 compares C m q' + Cmi compu¬ 
tations with data for a 9-caliber length body. GE-Spinner refers to the empirical 
curve fit of Reference 20. 

The final body computational example is taken from Reference 15. xk and xd 
are cone lengths, L, for an elliptic cone. Semi-minor to semi-major axis ratio is a/b = 
.75. a/L = .0866, x’ = 0. Figure 8 compares C\q' + C\i first- and second-order 
computations with that of Reference 15. Figure 9 compares Cmq' + C m i. Compu¬ 
tations for C.\or show significant differences between first- and second-order order. 
The computation of Reference 15 is a first-order computation using different numer¬ 
ical methods. Three computational planes seems to be adequate for this case. 


4.2 BODY-TAIL CONFIGURATIONS 

The first configuration considered for this section is the Basic Finner of Figure 
10 . 21,22 Figure 11 shows a comparison of data with computation for the roll damping 
derivative. The loading for the roll damping problem is low near the body and 
increases with span distance. The roll moment loading increases even more rapidly 
with span distance. At a side edge, the loading drops to zero. Load integration 
routines were modified to try to account for the sharp drop off of loading near a side 
edge. This improved C| p ’ prediction at lower Mach numbers. Figure 12 depicts a 
Magnus moment derivative, C.\ p -r, computation comparison with data. A body- 
alone estimate using the empirical computation of Reference 20 indicates a variation 
of-.4 to 3 for the Mach number range of 1.2 to 3. The computation here only follows 
for the data trend but not the magnitude. Figure 13 depicts a total pitch damping 
comparison. A trend of being above the data at lower Mach numbers is indicated 
here. 
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The next configuration is that of Figure 14(b) taken from Reference 23. The 
equation of the body in calibers is 

r b = -511 - (1 - x/6.25) 2 ! 3 ' 4 0 < x < 6.25 (4-2) 

= .5, x > 6.25 

Since drb/dx and 02 rb/ 0 x 2 are singular at x = 0, the nose is approximated as 


r= ‘29716844x - 131673491.x' 2 + 034340881x 3 0 <X<1 (4-3) 

b 

Equation (4-3) matches zero through second derivatives of Equation (4-2) at x = 1. 
The total pitch damping comparison is shown in Figure 15. Reference area is the 
extended to centerline fin planform area (two fins). Reference xr = 2xj) = the mean 
aerodynamic chord, 


C' = 


2 

3 


<C 2 + C 2 + CC) 

r t t r 

c + c 

l r 


(4-4) 


Here, C r and Ct are root and tip chords of the extended fin. The moment center is .2C' 
from the apex of the wing extended to the body centerline. Leading and trailing 
edges are supersonic and the computation compares fairly well with the data. 
Computations were adequate for the configuration of Figure 14(a) for supersonic 
leading edges, but very poor for subsonic edges when the moment center is .35C' from 
the extended wing apex. The reason is the inaccurate computation of C m i. For 
moment centers which are, more typically, not close to the point of action of C\i, the 
sensitivity and accuracy is better. 

The next configuration considered is depicted in Figure 16. A Cmq’ + C ma 
comparison for this configuration is shown in Figure 17. Here, the trend of over¬ 
prediction is reversed. 

The next body-tail configuration considered is the flechette of Figure 18 taken 
from Reference 24. Reference 24 (unpublished) compares one of the routines for 
computing total pitch damping in Reference 1 with available data for a large number 
of body-alone and body-tail configurations. The XM-144 body end radius of .005 
calibers is approximated as .05 calibers. Figure 19 compares computations with data. 
For this case, the leading edge of the fins is subsonic below Mach numbers of 3.25. 
The data was surprisingly smooth for this configuration. The boat tail has a 
significant effect on the computational results. 

The final body-tail configurations 0 f Figure 20 has six fins and a boat-tail 
angle of 2 deg. Figure 21 compares C| p ’ potential and local computations with the 
PNS computations of Reference 25. The sharp rise close to sonic leading edge 
conditions is predicted by the thin-fin, fin-alone methods of Reference 1. The boat tail 
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significantly affects the solution. Figure 22 shows a pitch damping computational 
comparison with data for a configuration close to that of Figure 20. The configuration 
and data are from Reference 24. Here, the conical nose length is 3.32 calibers and the 
moment center is 7.14 calibers from the nose. From Figures 21 and 22, one can see 
that the local solution is of limited value. 


4.3 BODY-WING-TAIL OR BODY-CANARD-TAIL CONFIGURATIONS 

Here, the additional interference of forward lifting surfaces on the tail is the 
primary phenomenon. As in earlier work, it is assumed that the velocity vector 
downstream of a supersonic trailing edge lies on a constant cylindrical ray plane, 0 
= 6f. 


The first computational comparison is for the B-C-T configuration of Figure 23. 
Figure 24 shows a computational comparison with experimental data. The trend at 
lower Mach numbers is as noted earlier. Here, the carryover canard to tail of C m i 
grows as the Mach number decreases. Evaluation of the plunging acceleration con¬ 
tribution of Equation (3-31) involves M^/p2 and differences which probably become 
sensitive as the Mach number is decreased. Computations of Ci p ’ also indicate a 
carryover effect which is not accounted for by component superposition methods. At 
Mo, = 1.76, C| p ' = -199.3 for the total configuration. For the same Mach number, the 
total of individual canard and tail configurations adds up to Ci p = -171.1. 

Figure 25 compares total pitch damping with data and a semi-empirical esti¬ 
mate for the RFL 122 configuration of Reference 26. The ordinate scale was not given 
and is inferred here. The C m « carryover phenomenon seems to be less extreme in this 
case. The semi-empirical estimate seems to be worse, assuming the data fit is correct. 
A Ci p ’ computation at M<r, = 1.5 yields -89.5 for the total configuration. The wing 
contribution to C| p at the same Mach number is -94.0. 

The final computation is for the B-W-T configuration of Reference 27. Total 
pitch damping comparison and configuration are depicted in Figure 26. Geometry 
and the moment center location are given in Reference 27. The geometry was scaled 
from Reference 27. x' is assumed to be at x' = 8.5 calibers from the nose. The C m i 
carryover phenomenon is noted as before. Missile refers to an Aerospatiale semi- 
empirical code. 
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5.0CONCLUDING REMARKS 


The computational methods, developed for computing static aerodynamic 
coefficients for noncircular bodies at supersonic Mach numbers, have been extended 
to the computation of dynamic derivatives. 

Geometric limitations are the same as for the earlier work. However, it is 
assumed that the bodies are pointed. 

The weakest elements in the original code are exaggerated here since moments 
are almost exclusively computed. The variation of pressures close to subsonic lead¬ 
ing and side edges is affected by edge singularities. C| p ' prediction is poorest for fins 
with low aspect ratios. C m q' + C m i prediction is poorest when the leading edge is 
subsonic and the point of action of C\ c .' + C\i is close to the moment center. The 
carryover effect for C m « is overpredicted when the leading edges are subsonic for 
bodies with two sets of lifting surfaces. 

The methods developed are capable of computing dynamic derivatives that are 
usually not computed or measured for in-plane and out-of-plane cases. 

Computational comparison with data is fairly good for most cases. In general, 
the repeatibility of dynamic derivative experimental values is much poorer than for 
static derivatives. 

The current computer code, implementing the dynamic derivative models, re¬ 
quires about 330,000 octal storage locations for a 15-by-60 grid for full plane 
computations. Computational times for axial rate computations are about the same 
as for static computations. Plunging rate computations are about 50 percent longer. 
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FIGURE 13. TOTAL PITCH DAMPING COMPARISON 
FOR THE BASIC FINNER 




(a) A • 3 SWEPT WING (3% BICONVEX CIRCULAR-ARC SECTIONS) A “ 3 UNSWEPT WING (3% BICONVEX CIRCULAR-ARC SECTIONS) 

FIGURE 14. ASPECT RATIO = 3 CONFIGURATIONS 
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FIGURE 15 PITCH DAMPING COMPARISON FOR AN 
AR = 3 WING-BODY CONFIGURATION 



OVERALL DIMENSIONS OF THE AF ATI AERODYNAMIC 
DATA CORRELATION MODEL 



FIGURE 16 AIR SHEW DEMONSTRATOR VEHICLE 
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FIGURE 21. ROLL DAMPING COMPARISON COMPUTATIONS 
FOR THE M735 PROJECTILE 



FIGURE 22 TOTAL PITCH DAMPING COMPARISON 
FOR THE M735 PROJECTILE 
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FIGURE 23. SIDEWINDER GEOMETRY 



FIGURE 24. TOTAL PITCH DAMPING COMPARISON FOR THE SIDEWINDER 
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FIGURE 25 TOTAL PITCH DAMPING COMPARISON FOR THE RFL 122 



FIGURE 26 TOTAL PITCH DAMPING COMPARISON FOR AN 
AEROSPATIALE MISSILE CONFIGURATION 
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NOMENCLATURE 

a Dimensionless speed of sound 

a B Reference speed of sound = free-stream value 

Ci General force or moment coefficient 

C p Pressure coefficient 

i Grid index for a constant x plane 

i' Unit vector in x direction 

j Grid index for constant B plane 

j' Unit vector in y or r direction 

k Grid index for constants plane 

k' Unit vector in z or 0 direction 

LM Aerodynamic symmetry mode 

Mi Moment or force 

M^ Free-stream Mach number 

n Unit normal vector from a solid surface 

p' Rotational rate about the x axis 

(NOTE: All q velocities are nondimensional) 
q Perturbation velocity vector 

q' Rotation rate about y axis 

q a Velocity vector associated with <{) a potential 

q. ix Velocity vector at t = 0 due to plunging acceleration 
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Velocity vector associated with 4>b potential 

Equivalent free-stream velocity vector relative to body axes 

Q - V, relative velocity 

Equivalent free-stream velocity vector relative to body coordinates at t 
= 0 

Free-stream acceleration relative to body axes (dimensionless) 

Equals q 0QO - i' ^ 

First-order velocity vector, steady part 
Total first-order velocity vector 
Second-order velocity vector 
First-order crossflow vector 

Total velocity vector of fluid in body axis coordinates relative to gas at 
rest 

1/2 p B V 2 , dynamic pressure 

Total first- or second-order velocity vector, t = 0 

Radial coordinate 

Rotation rate for z direction 

Body radius 

Body radius to outer fin edge 

Reference area 

Time 

Fin slope in x direction 
x component of perturbation velocity 
x component of c. g. acceleration 
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u', x component of q^ 

v r component of perturbation velocity 

v 0 y component of c. g. acceleration 

V Equals -q 0r ,, velocity relative to gas at rest in body axis coordinates 

Vuo Magnitude of velocity vector of c. g. at t = 0 

w 0 component of perturbation velocity 

w 0 z component of e. g. acceleration 

x Cartesian coordinate from nose — zero incidence direction 

x' Moment center distance from nose 

X|{ Reference length 

xj) Reference length associated with reduced frequency 

y Cartesian coordinate perpendicular to vertical body symmetry plane (to 

the right looking downstream) 

z Cartesian coordinate perpendicular to x, y plane (up, looking down¬ 

stream) 

« Angle of attack 

P Equals VmJ- 1 

P’ Angle of side slip 

Y Specific heat ratio 

£q Equals (dr^/fjO) / r^ 

C, Equals ^ (£,), clustering transformation in £, direction 

0 Cylindrical coordinate angle from leeside plane 

E, Equals (r - r^Vlx/p - r t) ), shearing transformation 

p Density 


A-5 








NSWC TR 89-99 


p r Free-stream density 

4> Disturbance potential 

(J>, First-order disturbance potential for q, = 0 

4>; Equals + (Moo/p 2 )^ + V^ft - xM 2 / (p 2 V 00 )]<|)t > 

(Total first-order disturbance potential $ a and 4>b are equivalent static 
problem potentials.) 

(f> 2 Second-order disturbance potential 

d> Equals (B) - clustering transformation in B direction 

X Absolute temperature ratio 

ojj General term in computing reduced frequency 
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